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We explore a nonlinear field model to describe the interplay between the ability 
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q ' of excitons to be Bose-condensed and their interaction with other modes of a crystal. 

We apply our consideration to the long-living para-excitons in C112O. Taking into 

account the exciton-phonon interaction and introducing a coherent phonon part of 

the moving condensate, we derive the dynamic equations for the exciton-phonon 

condensate. These equations can support localized solutions, and we discuss the 

conditions for the moving inhomogeneous condensate to appear in the crystal. We 

calculate the condensate wave function and energy, and a collective excitation spec- 

q^ ■ trum in the semiclassical approximation; the inside-excitations were found to follow 

ON . the asymptotic behavior of the macroscopic wave function exactly. The stability 

conditions of the moving condensate are analyzed by use of Landau arguments, and 

Landau critical parameters appear in the theory. Finally, we apply our model to 

describe the recently observed interference and strong nonlinear interaction between 

two coherent exciton-phonon packets in CU2O. 
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PACS numbers: 71.35.H-z, 71.35.Lk 



1 Introduction 

Excitons in semiconductor crystals [||] and nanostructures @j,|| are a very interesting 
and challenging object to search for the process of Bose-Einstein condensation (BEC). 
Nowadays there is a lot of experimental evidence that the optically inactive para-excitons 
in CU2O can form a highly correlated state, or the excitonic Bose Einstein condensate 
0,0,0- A moving condensate of para-excitons in a 3D CU2O crystal turns out to 
be spatially inhomogeneous in the direction of motion, and the registered velocities of 
coherent exciton packets turn out to be always less, but approximately equal to the 
longitudinal sound speed of the crystal B. 



Analyzing recent experimental H,@,[0] and theoretical |8|]-[T3[ studies of BEC of 
excitons in CU2O, we can conclude that there are essentially two different stages of this 
process. The first stage is the kinetic one, with the characteristic time scale of 10 ~ 20 ns. 
At this stage, a condensate of long-living para-excitons begins to be formed from a quasi- 
equilibrium degenerate state of excitons (// 7^ 0, T c s > T\ att ) when the concentration 
and the effective temperature of excitons in a cloud meet the conditions of Bose-Einstein 
Condensation 0. Note that we do not discuss here the behavior of ortho-excitons (with 
the lifetime r ort ho — 30 ns) and their influence on the para-exciton condensation process. 
For more details about the ortho-excitons in CU2O, ortho-para-exciton conversion, etc. 
see f!J 



14 , 15 



The most intriguing feature of the kinetic stage is that formation of the para-exciton 
condensate and the process of momentum transfer to the para-exciton cloud are happening 
simultaneously. If the diameter of an excitation spot on the crystal surface is large enough, 
S'spot — S'surf, and the energy of a laser beam satisfies e p hot 3> E gap , nonequilibrium 
acoustic phonons may play the key role in the process of momentum transfer. As a result, 
the mode with macroscopical occupancy of the excitons appears to be with (k) 7^ 0, where 
h(k x ) = m x v and v is the packet velocity. 

Indeed, the theoretical results obtained in the framework of the "phonon wind" model 
I(|,[|l6j and the experimental observations @],||,|| are the strong arguments in favor 



of this idea. To the authors' knowledge, there are no realistic theoretical models of the 
kinetic stage of para-exciton condensate formation where quantum degeneracy of the 
appearing exciton state and possible coherence of nonequilibrium phonons pushing the 
excitons would be taken into account. Indeed, the condensate formation and many other 
processes involving it are essentially nonliner ones. Therefore, the condensate, or, better, 
the macroscopically occupied mode, can be different from n(k = 0) ^> 1, and the language 
of the states in k-space and their occupation numbers n(k) may be not relevant to the 
problem, see [fL7| . 

In this study, we will not explore the stage of condensate formation. Instead, we 
investigate the second, quasi-equilibrium stage, in which the condensate has already been 
formed and it moves through a crystal with some constant velocity and characteristic 
shape of the density profile. In theory, the time scale of this "transport" stage, At tr , 
could be determined by the para-exciton lifetime (r para ~ 13 /xs [Q). In practice, it 
is determined by the characteristic size £ of a high-quality single crystal available for 



experiments: 

A£ tr - tjc x ~ 0.5 - 2 fis < r para , 

where q is the longitudinal sound velocity. 

We assume that at the "transport" stage, the temperature of the moving packet (con- 
densed + noncondensed particles) is approximately equal to the lattice temperature, 

T e S = ^latt < T c . 

Then we can consider the simplest case of T = and disregard the influence of all sorts of 
nonequilibrium phonons (which appear at the stages of exciton formation, thermalization 
|10| ) on the formed moving condensate. 

Any theory of the exciton BEC in CU2O has to point out some physical mechanism(s) 
by means of which the key experimental facts can be explained. (For example, the con- 
densate moves without friction within a narrow interval of velocities localized near q, and 
the shape of the stable macroscopic wave function of excitons resembles soliton profiles 
0.) Here we explore a simple model of the ballistic exciton-phonon condensate. In this 
case, the general structure of the Hamiltonian of the moving exciton packet and the lattice 
phonons is the following: 

H = H cx {ft, •$) - vP ex (^, j>) + H ph (u, ft) - vP ph (u, ft) + flintftW, djU k ). (1) 

Here i\) is the Bose-field operator describing the excitons, u is the field operator of lattice 
displacements, ft is the momentum density operator canonically conjugate to u, and P is 
the momentum operator. Note that the Hamiltonian (|I|) is written in the reference frame 
moving with the exciton packet, i.e. x — > x — vt and v = const is the ballistic velocity of 
the packet. 

2 3D Model of Moving Exciton-Phonon Condensate 

To derive the equations of motion of the field operators (and generalize these equations 
to the case of T 7^ 0), it is more convenient to start from the Lagrangian. In the proposed 
model, the Lagrangian density has the following form in the co-moving frame 

£ = -J ($ dt $ ~ dt $ W + v ~n ( dx * ^ ~ ^ dx ^~ 

-4W - |WW> - |(^(x,t)) 2 (^(x,t)) 2 - |(^(x,t)) 3 ^(x,t)) 3 + 

+f (d t u) 2 - ^-{d^f - P A Kz{dj u s f -^(d t ud x u + d x iid t u) + ^(d x h) 2 - 

-^ t (x,^(x,t)Vu(x,{), (2) 

where m is the exciton "bare" mass (m = m e + rrih ~ 3m e for Is excitons in CU2O), uq is 
the exciton-exciton interaction constant (uq > corresponds to the repulsive interaction 



between para-excitons in CU2O p!S[), p is the crystal density, <To is the exciton- longitudinal 
phonon coupling constant, and v = (v, 0, 0). The energy of a free exciton is E g + h 2 'k 2 /2m. 
Although the validity of the condition n a| <C 1 (5b is the exciton Bohr radius) makes 
it possible to disregard all the multiple-particle interactions with more than two partici- 



pating particles in H cx [19|, we include the hard-core repulsion term originated from the 



three-particle interaction in C, i.e. v\ 7^ and 



< v\ja\ <C ^o/ a B — con st Ry* 



(For 3D case, one has to take const ~ 10 because u = An (h 2 /m) a sc and a sc ~ (1 ~ 3) a-Q] 
see the discussion in |pCfl .) 

Moreover, in the Lagrangian of the displacement field, we include the first nonlinear 
term oc K3 (du) 3 . (The dimensionless parameter K3 originates from Taylor's expansion of 
an interparticle potential U(\ri — rj\) of the medium atoms.) Assuming that a dilute 
excitonic packet moves in a weekly nonlinear medium, we will not take into account more 
higher nonlinear terms in ([|). 

For simplicity's sake, we take all the interaction terms in £ in the local form and 
disregard the interaction between the excitons and transverse phonons of the crystal. 
Note that the ballistic velocity v is one of the parameters of the theory, and we will not 
take into account the excitonic normal component and velocity, i.e. v = v s ~ Vv? c > 
(T = 0). This means that we choose the spatial part of the coherent phase of the packet, 
(^ c (x), to be in the simplest form, 

exp(i<£> c (x)) = exp(i(<£> + k x)), ip = const, hk = mv. (3) 

The equations of motion can be easily derived by the standard variational method 
from the following condition: 

SS = 6 f dtdx.C($(x.,t), ^(x,£), u(x,t)) = 0. 

Indeed, after transforming the Bose-fields ip' and ip by 

"0(x, t) — ► exp(—iE g t/h) exp(imvx/h)ip(~x.,t), 
we can write these equations as follows: 

(ihd t + mv 2 /2)^(x,t) = 

h 2 

= ( - —A + z/ V^(x, t) + z/i^ t2 ^ 2 (x, t) )^(x, t) + a Vu(x, t) $(x, t), (4) 

(d 2 - c 2 A - v(d t d x + d x 8 t ) + v 2 d 2 x )u s (x, t) - c 2 5> 3 d 2 u s a,-« a (x, t) = p~ V d s (^(x, t)). 

3 

(5) 

We assume that the condensate of excitons exists. This means that the following 

representation of the exciton Bose-field holds: ip = tpo + 8ip. Here tpQ 7^ is the classical 



part of the field operator ip or, in other words, the condensate wave function, and 8ip is 
the ffuctuational part of if;, which describes out-of-condensate particles. 

One of the important objects in the theory of BEC is the correlation functions of 
Bose-fields. The standard way to calculate them in this model (the excitonic function 
(■?/>(x, t)^'t'(x / , £')), for example,) can be based on the effective action or the effective Hamil- 
tonian approaches |^TJ. Indeed, one can, first, integrate over the phonon variables u, get 



the expression for S e ff (ip, ijr) and, second, use S e s (or H cfi ) to derive the equations of 
motion for i[) , Si[>, correlation functions, etc.. 

In this work we do not follow that way; instead, we treat excitons and phonons equally 



TTJI , [£2| , |23[ . This means that the displacement field u can have a nontrivial coherent part 
too, i.e. u = uo + Su and Uo 7^ 0, and the actual moving condensate can be an exciton- 
phonon one, i.e. ^(x, £) ' u o( x , £). Then the equation of motion for the classical parts of 
the fields ip and u can be derived by use of the variational method from C = C(ip, ip*,u), 
in which all the fields can be considered as the classical ones. Eventually we have 

(ihd t + mv 2 /2)ip (x,t) = 

h 2 

- — A + i/ |^o| 2 (x,t) + za|^ | 4 (x,£) Vo(x,£) + cr Vu (x,t)^o(x,t), (6) 

2m J 

(d 2 -c 2 A-2vd t d x +v 2 dl)u 0s (x,t)-c 2 J2 2K 3d]uo s d j u 0s (x,t) = p _ V 9 s (|^ | 2 (x,t)). (7) 

j 

Notice that deriving these equations we disregarded the interaction between the classical 
(condensate) and the ffuctuational (noncondensate) parts of the fields. That is certainly 
a good approximation for T = and T <^iT c cases |]24}| . 

In this article, a steady-state of the condensate is the object of the main interest. 
In the co-moving frame of reference, the condensate steady-state is just the stationary 
solution of Eqs. (|6|), (0) and it can be taken in the form 

V>o(x, t) = exp(-i/it) exp(iy3)0 o (x), u (x, t) = q D (x), 

where <p and q are the real number functions, and if = const is the coherent phase of 
the condensate wave function in the co- moving frame, see Eq. (^). (This phase can be 
taken equal zero if only a single condensate is the subject of interest.) 
Then, the following equations have to be solved, (fi = jl — mv 2 /2 ): 

h 2 

M (x) =( - —A + z/ o 0o(x) + ^i0o(x) )0 o (x) + (ToVq„(x) o (x) , (8) 

- { (c 2 - v 2 )d 2 x + c 2 d 2 + c 2 d 2 }g os (x) - c 2 £ 2« 3 d 2 q os ^g os (x) = p-V o ^0 2 (x). (9) 

Note that in order to simplify Eq. (0) to Eq. (|9D, we assumed only m (x, t) = q (x). In 
this model, it is enough to obtain localized solutions for the displacement field. 



3 Effective ID Model for the Condensate Wave Func- 
tion 

Solving Eqs. (§),(0) in 3D space seems to be a difficult problem. However, these equations 
can be essentially simplified if we assume that the condensate is inhomogeneous along the 
x-axis only, that is 

o (x) = 4> {x) and q (x) = (q {x), 0,0). 

Note that the cross-section area of an excitation spot, S, has to be basically constant across 
the sample cross-section. In this case, the problem can be considered as an effectively 
one-dimensional one. 

Such an effective reduction of dimensionality transforms difficult (nonlocal differential) 
equations for the condensate wave function into a rather simple differential ones, and 
obtained in this way the effective ID model for the condensate wave function O • q Q 
conserves all the important properties of the "parent" 3D model. 

Indeed, if v < q, the following equations stand for the condensate (y(x) = d x q (x)): 

p,<f> (x) = ( -{h 2 /2m)d 2 x + z/ o 0o(x) + za0oO) )0oO) + o"o y(x) (f>o(x), (10) 

- (c 2 - v 2 )d x y{x) - 2c 2 n 3 d x yy(x) = p" 1 a d x <j> 2 {x) . (11) 

The last equation can be easily integrated, 

y(x) + k 3 y 2 (x) = $(x) + const, (12) 

and solved relative to y(x). Here 

2 
^3 = 2 ' 2 «3 = l{v) «3, $(X) = 2 - <f) 2 (x) = -7(v) — ^ </>l(x). 

cf-v 2 p{cf-v 2 ) pcf 

Note that the medium nonlinearity parameter k 3 can be enhanced by the factor of the 
order of 4 ~ 10 if the value of v is less, but close to q. (For spatially localized solutions, 
d x qo{%) — and (f> 2 (x) ~ at \x\ ^> L c h, so that const = 0.) 

If «3 < 0, we can always represent the solution of ([Til ) , ([12]) i n the following form 



y (x) = $( x ) + |k 3 | $ 2 (x) + 2 k\ $ 3 (x) + • ■ • . (13) 

(Indeed, the parameters of medium nonlinearity can be chosen as k 3 < and k 4 > 
5[.) After substitution of (p~3|) into Eq. (|T0"D, Eqs. (|10|), ([TT]) can be rewritten as follows: 



/i0 o (x) = (-(/ i 2 /2m)^ + ^^( 2; )+^^(x) + e 2 )0 o (x), (14) 

d x q Q (x) = ®(x) + \R 3 \<S> 2 (x) + e' 2 , (15) 

where the interparticle interaction constants are renormalized as follows 

u = Vq - a —- -, vx = v x + a -—— — -j « 3 , (16) 

p{cf-V 2 ) p (^_^2 



and more higher nonlinear terms are designated by e-i- A small parameter in Eq. ( IB ) 
comes from the term 

a /p(c 2 -v 2 ) = -y(v)(a /Mct)al 

where 7(f) = c 2 /(c 2 — v 2 ) and M is the mass of the crystal elementary cell. 

The effective two-particle interaction constant vq (v) can be negative if the velocity 
of the condensate lies inside the interval v < v < q , where 



^o = Jcf - (o-g/pi/o). (17) 



Outside this interval, vq (v) > JT| and the velocity v can be called the first 'critical' 



velocity in the model. The meaning of this velocity can be clarified by rewriting fll6|) in 
the dimensionless form, 



lip) TT5 , (1* 



If v > v, 



O' 



cto \ ^0 const Ry* a 3 

Mtf) > a af ~ o~ ~ 3 ' 



7(«) 773 > 7-3 * —^ 77 > (19) 



where Ry* and a 2 B are the characteristic energy and the cross-section of two-particle 
collisions in the exciton subsystem. The following inequalities are true for excitons in a 
crystal 

a|>(>)a 3 and const Ry* < (<) <J , 

and, usually, the value of the parameter vqJgq af is > 1. 

For para-excitons in Cu 2 0, however, we assume the (effective) value of vq/o-q af can be 
estimated as 0.3 ~ 0.6 < 1, whereas the value of ao/Mc 2 ~ 0.1 ~ 0.3. This makes the 
inequality (|19D valid at, say, v ~ (0.8 ~ 0.9) q, or 7(f) ~ 5. Thus, within the effective 
ID model, the critical factor 7 = 7(1*0) is the following ratio 

70 = y^aj) ' \t£%) ' 

and, for the substances with uo/aoaf < 1, the regime with uq < can be obtained at 
velocities reasonably close but not equal to q, for example, beginning from some 7 < 10, 
( 7 (0.95q)~10). 

On the other hand, the effective three-particle interaction constant V\ (v) is always 
positive for crystals with K3 < 0. It can be represented in the dimensionless form as 
follows 

^33 - -nst' f^O % + j(v) |«3| Lv) *)'. (20) 

ffoK) 2 \o-oafj af V Mc i) 

Here we estimated the "bare" vertex of the three-particle collisions as 

V\ ~ const Ry*a^ ~ const' uq a B , const' < 1, 
and the same Ry* can be taken as a characteristic energy of collisions. The effective vertex 



v\ > is enhanced by the medium nonlinearity, and the both terms in the r.h.s. of ([20]) 
can be equally important at 7(f) > -y Q . 



Note that in the case of strongly nonlinear lattices with excitons, the effective interac- 
tion vertices in ( |14j ) ( u{, £% etc.) depend on the velocity v and the parameters of medium 
nonlinearity (k^, Ka, etc.). Then the effective exciton-exciton interaction can be strongly 
renormalized at sufficiently large gamma- factors 7(f) and the vertices may change their 
signs as it can happen with Vq (v). In this article, however, we consider the case of weakly 
nonlinear medium with excitons (e.g., a crystal with long living excitons). More accu- 
rately, this means that at velocities v — > c\ the effective vertex v§ (v) became < 0, while 
the more higher vertices, such as 17 (v) and v^. (v), do not change their sign; they remain 
> at j(v) > 7o- Finally, to describe the weakly nonlinear case, it is enough to take into 
account the parameters V\ > and K3 < and neglect more higher nonlinearities (e2 
and e' 2 in flU, ©). 

In this study, we will consider the case of v < v < Ci in detail. Indeed, in the case 
of Vq{v) < and 17(f) > 0, some localized solutions of Eqs. (|14|) , (|T5|) do exist. For 
example, the so-called 'bright soliton' solution of ([14]) exists if the generalized chemical 
potential is negative, fx < 0, and \p\ < //*. Here 

P = C\r/q\~ ~0.2<7 — — g-r-. (21) 

(16/3) Z7 (vi/(Toa?) 

For |k 3 | ~ 1 and 7(f) > 7 ~ 3 ~ 5, we can roughly estimate the effective vertex v\ (v) 
as 

^(u)/<t o? ~ (1 ~ 10) {vo/<T af). 

Then n*(v) ~ (10 _1 ~ 1CT 2 ) Ry*, and the more is the value of | ^3 1 the less is the value 
of fi*(v). 

The 'bright soliton' solution of Eq. (]14]) can be represented in the following form 



&(s) = $o/(/?($o)z,m(*o)), W = \l— lAK^o). (22) 



'2m 

7? 



Here r/i($ ) is some dimensionless parameter, and the generalized chemical potential 
pi, < is given by the formula 

| / i| = | / i|($ ) = |^|^/2-^^/3. (23) 

Like the chemical potential |//|, the amplitude of the bright soliton, $ , satisfies 

$ 2 <TO 2 = l^l/(4/3^), 

and /i* = |i^| $* 2 /4. 

For |/2|//i* <C 1, the following approximation is valid 

»7i«^(l/i|/A«*) + ^(l/i|//) a <l, (24) 

and this formula can be used up to \ji\/ n* — 0.5. Then we can approximate the solution 
of (0) by the following formulas 



4o(x) « $0 ( Vl - r/i($ ) cosh(/5($ ) x) + (1 - ^1 - ^($0)) ) , (25) 



o (x)~2$ o exp(-/?($ o )|x|)/^l-r ?1 for \x\ > 2/?($ )~\ \(i\ < //. (26) 

The amplitudes of the exciton and phonon parts of the condensate, the characteristic 
width of the condensate, and the value of the effective chemical potential \x depend on the 
normalization of the exciton wave function <p (x). We normalize it in 3D space assuming 
that the characteristic width of the packet in the (y, z)-plane is sufficiently large, i.e. the 
cross-section area of the packet S± can be made equal to the cross-section area S of a 
laser beam and 

o_l — o — O sur f. 

Then we can write this condition as follows: 

yjVo| 2 (M)dx =SJ<j>l{x) dx = N , (27) 

where N Q is the number of condensed excitons, and, generally, N ^ N tot . 
Applying this normalization condition, we get the following results 

°^ 2{N*/N fxRj*~a% + 217 l {v) 

Here we used the following notations, N* = 2S/a%, h 2 /m = 2xRy*a 2 3 , where Ry* = 
h 2 /2[i exc a 2 3 and x = fi cxc /m. The formula fl2"8|) is valid for \fi\ffj,* < 0.3 ~ 0.4. We 



assume that, at N*/N Q = n Q > 10 (this is the important parameter!), we always have 

2n 2 xRy*a 6 B > u{{v) = ei(|« 3 H 4 - (1 ~ 10) Ry*a|- 
Then, the following inequalities are valid: <&|(JV , v) -C $* 2 and 







MNo, 


v) ~ 

; 2{2n 2 x 


H(v)\ 2 
Ry*a| + 4 


"i (v) } 


__ H(v)\ 2 
5.3 v\ (v) 






(29) 


** 


The characteristic 
(v) =e (v)a%): 


length of the 


packet can 


be estimated from Eq. 


m 


as 


follows 



^(Ar„, )s! !«•„)«! IWI : 



4 " 4(2xRy*a B ) ^2 + (i^j/^Ry*) 
1 |e (v)| 1 



p , _ — at n o >10. (30) 

8 x Ry a B n 

Therefore, at j(v) ~ 2j Q , we can roughly estimate 

Rv* 

L ch (iV , u) ~ 4 — -f- n a B ~ 4 (10 1 n ) o B , (31) 

lw)l 

and, for the average concentration of condensed excitons in the packet, n Q , we have 

n a% « (A^ /5L ch ) o| ~ 1/n 2 < 1. 

9 



Recall that the second part of the condensate, the displacement field q (x), is of the 
same importance as the first part, the exciton wave function 4> (x). The displacement 
field d x q (x) can be represented as follows 



2 



d x q (x) = -j(v) (-^J (af <Pl(x) ) + j(v) |«s| [l(v) -j^j (af f {x) ) . (32) 

To estimate its amplitude, d x q , we have to estimate the parameter af $ 2 first. 
For n > 10 and |e~ (v)\ ~ (10 _1 ~ 1) Ry* (i.e. 7(f) > 2^ ), we obtain 



af \e (v)\ 1 _ af 1 
a% xRy* 2n 2 a% 2n 2 Q 



u i ^o — ~3 „r>-.* r>-o ~3 0^2 ^ Jv o' 



If this parameter is small enough, such as af $ 2 (iV , v) ~ 10 3 ~ 10 5 , we can neglect 
the nonlinear corrections to the amplitude d x q Q < and to the shape of d x q (x) as well, 

d X 1o = -l(v) ^ (af $ 2 ) {l - 7 („) \k 3 \ ( 7 („) ^) (af $ 2 ) } « -7(«)(^)(a? ^)- 

(33) 
Thus, due to the validity of n a% -C 1, there is almost no difference between the approx- 
imation 

4>o{x) ■ d x q Q (x) « $ cosh- 1 (/5($ ) a;) ■ (-\d x q \ ) cosh- 2 (/3($ ) x), (34) 

where 

^(^^Lo- 1 ^^^— V, iV = iV 7n , n »l, (35) 

x Ry 2 n a^ 

and the exact solution of the weakly nonlinear case with V\ > and k 3 < 0. For S ~ 
(10~ 2 ~ 10~ 3 )cm 2 , a| ~ (25 ~ 50) 10~ 16 cm 2 , we estimate N* ~ 10 13 ~ 10 14 . Although 
the approximate solutions we used in this study are valid for N <C N*, they can be used 
at N Q < N* for estimates. 

Note that the effective chemical potential is a rather small parameter in this model, 

4 n* x Ry a^ 4 n^ \ x Ry / 



That is why the characteristic length, L c h oc |/i| _1 / 2 , see fl3"0|), can be estimated as (10 2 ~ 
10 4 ) as within the validity of approximations (|28|), (|29|). Moreover, |/i|/ju* < 10~ 2 and 
the parameter 771 (<& ) in (^) can be estimated as ~ 10~ 2 . In this case, one can neglect it 
in Eq. (p5|) . 

Returning to the laboratory reference frame, we can write the condensate wave func- 
tion in the form (see Fig. 1): 



mv 2 



ipo(x, t) ■ u (x, t)5\j ~ exp I —i I E g -\ |//| \ t \ exp(i(ip + mvx) ) x 



10 



x fcoCOslf^Lo^x-v*))- (Qo - QotanhiL^ix-vt))), (37) 



where we count the exciton energy from the bottom of the crystal valence band, (E g < 
i? ga p), and 2Q (N , v) is the amplitude of the phonon part of condensate, 



(^f)(ii) a,<< 



Qo ~ l{v) (-t^-o ) \-zk- — ) ai < a t 



To calculate the energy of the moving condensate within the Lagrangian approach, 
(see Eq. (||) ), we have to integrate the zeroth component of the energy- momentum tensor 
Tq over the spatial coordinates. Consequently, we have the following formula 

r °(x, t) = 4^> + ^ vc V0 O + 1 (0:) 2 </> 2 + 1 (^) 3 ^+ 

+ |(«9 4 g ) 2 + ^ (^.g ) 2 H- ^ « 3 (<9,g Q ) 3 + <x 0> o 0,g o . 

Here we do not take into account a small correction to this energy due to the quantum 
depletion of the condensate ((5\E^<5\I/(x)) T=0 7^ and ((d x 6uj) 2 )T=o ¥" 0)- Then the 
result reads 

E (N , v) = |rfxT ° = E cx + E- mt + E p 






,iV o ^H-^J-iV o (|A|H-^0^3)H-iV o |^_r^7>) ^J ~ (a? *> 

(38) 
We will disregard the terms ~ A^z/i^ in i? int < 0, and the corrections oc \k 3 \ in .E7 p h. 
Then we can write 

\E- mt \/N « |i£(i;)|$*/2 + z/o$ 2 /3 ~ (V4) (4$o) < %* ( 39 ) 

and 

M (r 2 4- v 2 ) 
E ph /N « U 2 + J tf(iV , «) ~ Mc? 0(JV o , v), (40) 

where 

Wo,v)=(>y(v)£ s ) |(af^)«l. (41) 



Note that the parameter $(N Q , v) is a rather small one, •& ~ af $ 2 ~ 10~ 3 ~ 10 -5 , so 
that the value of E ph /N can be < Ry*, and, roughly, $ 2 oc N 2 . 

One can see that the exciton-phonon condensate carries a non-zeroth momentum, 

mix -'ex, a; "T -*ph, a;- 

P oa; = dx(h/2i)((j)* (x,t)d x (j)o(x } t) -d x (f)* (x,t)(f)o(x,t)) - pd t u (x,t) d x u (x,t) = 

dxmv 4> 2 (x) + ^(tC^tft a iVo( x ) ) ~ ^mf + N Mvd(N , v) = 

= JV m{l + (M/m)${N , v) } v. (42) 

Thus, we obtain m c g = m {1 + (M/m)d(N , v) } and estimate the parameter 
(M/m) $(N , v) ~ 1 ~ 5 at 7(1;) > 27 , n > 10. 
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4 Low-Lying Excitations of Exciton-Phonon Conden- 
sate 

To consider the stability of the exciton-phonon condensate moving in a lattice, one has 
to couple the excitons with different sources of perturbation, such as impurities, thermal 
lattice phonons, surfaces, etc.. In this work, however, we will not specify any source. 
Instead, we consider the stability conditions in relation to creation (emission) of the 
condensate excitations that can be found in the framework of investigation of the low- 
energy excitations of the condensate itself. 

Although the condensate wave function (j) {x) ■ q (x) was obtained in the framework of 
the effective ID model, we normalized it in 3D space. Therefore, we can use this solution 
as a classical part in the following decomposition of 3D field operators in the co-moving 
frame: 

■0(x, t) = exp(-i/it) (0 o (x) + £0(x, t) ), (43) 

Uj (x, t) = q Q (x) Sij + 5uj (x, t) , (44) 

where fi = jl — mv 2 /2. Substituting the field operators of the form (f43"D , (|44l) into the 
Lagrangian density (H), we can write the later in the following form: 

£ = £ o (e"^0 o (x), q (x)6 lj ) + £ 2 (#t(x,t), #(x,t), 6u(x,t)) + ..., (45) 

where C stands for the classical part of C, and Li is the bilinear form in the ^-operators. 
In the simplest (Bogoliubov) approximation [Effl,p7|, £ « C Q + C 2 and, hence, the 



bilinear form £ 2 defines the equations of motion for the fluctuating parts of the field 
operators. As a result, these equations are linear and can be written as follows: 

ihd t <V(x, t) = I -g^A + |/2| + {is + Mv)} 4>l(x) + {2v x + V^v)} <j>i(x) J 6$(x, t) + 

+ (uo<f>l(x) + 2v 1 <ft(x) ) 5^(x, t) + a o 0o(x)V5u(x, t), (46) 

and 

(d 2 - cfA - v(d t d x + d x d t ) + v 2 d 2 x ) %(x, t) = 

= p-V o ^-(0 o (x)(^(x,t) + ^t( x ,t))), j = 2,3(=L), (47) 

(d 2 t - c 2 A - v(d t d x + d x d t ) + v 2 d 2 x ) 6u x (x, t) - 
-c 2 2k 3 (d x q Q (x)) d%6u x (x, t) - c 2 2k 3 {d 2 x q {x)) d x 5u x (x, t) = 
= p"V o ^.(0 o (x)(^(x,t) + ^ t (x,t))), j = l(=x). (48) 

The same approximation can be performed within the Hamiltonian approach. Indeed, 
decomposition of the field operators near their nontrivial classical parts leads to the de- 
composition of the Hamiltonian ([!]) itself, and - as it was done with the Lagrangian - 
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u (x-vt) 



(|) (x-vt) 




exc & phon 



x 



■> 



1-inside-exc 



Figure 1: Moving exciton-phonon condensate, as it appears in the quasi-stationary 
model, (p (x — vt) • u {x — vt)Sij, is presented by bold lines on this Figure. Here, 2Q Q is 
the amplitude of the coherent phonon state u (x — vt), and $ is the amplitude of the 
macroscopic wave function of excitons. Longitudinal exciton-phonon excitations (k 1 1 Ox) 
of the condensate are schematically depicted. Under transformation iV — » N — 5N, the 
condensate wave function is changed as it is presented by dashed lines. 
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only the classical part of H, H , and the bilinear form in the fluctuating fields, H2, are 
left for examination: 

H « H o (^* , Vo, vr , M ) + H 2 (5ft, 8$, 8w, 8u). (49) 

In the co-moving frame, fij = pd t Uj — pvd x itj, i.e. 

7r j(x) = — pv d x q Q (x) Sij 7^ and Sitj = pd t 5iij — pv d x Suj, 

and the standard commutation relation, [5uj(x., t), <$7T a (x',£) ], has the form 

[5wj(x, t), pd t Su s (x.',t) — pvd x Su s ('x',t)} = ih5('x — x.') 5j S . (50) 



However, the Hamiltonian (|49[) can be diagonalized and rewritten in the form: 



H = H Q (e %fJ,t (f> (x), q (x) ) + 5E Q + ^ hu 1>s &! )a ai )3 + X! ^M al^cb.s- (51) 



1,8 2,s 



Here, 5.E is the quantum correction to the energy of the condensate and the indexes 1, s 
and 2, s label the elementary excitations of the system. We assume the operators d] >s , 
d JiS are the Bose ones. These operators describe two different branches of the excitations, 
j — 1, 2, and they can be represented by the following linear combinations of the "delta- 
operators" : 

& hS = |dx ( U j)S (x) <ty>(x) + VJ- s (x) <tyt(x) + Y* t8 (x) 6ut(x) + 4» 5# *( x ) ) ' ( 52 ) 
4, = |dx ( £/*» # f (x) + V£(x) #(x) + y£(x) <^(x) + Zj» <^(x) ) . (53) 



Note that by analogy with the exciton-polariton modes in semiconductors |28| , |29[ the 
excitations of the condensate ( |3"T| ) can be considered as a mixture of exciton- and phonon- 
type modes. However, in this model, the phonons are fluctuations of the (n (x, t), Uo(x, t) ) ■ 
part of the condensate. The commutation relations between a-operators are Bose ones, 
so that 

lead to the following orthogonality condition 

fdx (u^uiA*) - v hs v* s ,(x) ) + m E A x (nX'( x ) - KsKAx) ) = **■ 

J r=l,2,3 J 

Since the a-operators (see Eq. (pT]) ) evolve in time as simply as 



a jtS (t) = Q~^ st a jtS , a} s (t) = e^- si a ] 



3,si 



these operators (and the frequencies {u>j iS }) are the eigenvectors (and, correspondingly, 
the eigenvalues) of the equations of motion (|46| ) ,(|4~7| ) obtained within the Lagrangian 
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method. Then, the time dependent "^-operators" in (|4"6j) , (|£7|) can be represented by the 
following linear combinations of the a-operators: 



6#{x,t) =^u M (x)« M e- i ^ i + v^(x)d} >8 e <Wl .' t + 

1,8 

+ E u 2, s (x)d 2 , s e-^-* + v^x)^^-', 
<5£ r (x, t) = 2 C[ )S (x) di, a e-^-* + C£(x) a^e*" 1 "' + 



(54) 



l,s 



^C 2 r s (x)d 2 , s e 

2,s 



-*■*-* + C2>) aLe**.-*, 



(55) 



For <57T r (x, £), one has to change CJ >s (x) to £)£ = p(—iujj iS — vd x )Cj S (x) in (|55"D. Note 
that this ansatz is, in fact, a generalization of the u-v Bogoliubov transformation. 
Then we can rewrite Eqs. fl52|) , (|53|) as follows (j = 1, 2) 

a j>s = Jdx (u*» #(x) - v* s (x) # f (x) - (*/n)£ft(x) 5n r (x) + {i/h)C£{x) <Wh r (x) ) , 

(56) 

4, = |dx (u iiS (x) #t( x ) - Vj> (x) <ty(x) + (i/ft)Z£.(x) ^ r ( X ) - (i/ft)q.(x) 57T r (x) ) , 

(57) 
and one of the orthogonality relations has the form (s = s') 

|rfx(|n liS (x)| 2 -|v M (x)| 2 ) + 

+ {i/h) J2 [dx(C[* sP (-iuj 1>s -vd x )Cl s (x) + p{-iu 1>s + vd x )C[* s Cl s (x)) = l. (58) 

The question we want to clarify is whether coupling between excitonic excitations and 
phonon excitations is important for understanding the condensate excitations. Substi- 



tuting ansatz (J54])-(|55D into Eqs. (f46|) , (j47|) , we obtain the following coupled eigenvalue 



equations [ |TT|j : 

(L(A) - hu j>s ) u,, s (x) + (mI(x) + 2M 4 Q (x)) v,, s (x) + a o 0o(x)a r q s (x) = 0, (59) 

(u 4> 2 o (x) + 2v l( f> 4 {x)) u ilS (x) + {1(A) + huo hS ) v i)S (x) + cr o o (a;)9 r q s (x) = 0, (60) 



-p i o- o <9 r (0 o (a;)u jiS (x)) -p V o a r (0 o (x)v iiS (x)) + 



+ 



-ia^ s - v d x ) 2 - c 2 A 



CL(x) = 0, r = 2,3, 



(61) 



+ 



-P V o 9 z (0 o (x)u J> (x)) -p V o 9 z (0 o (x)v iiS (x)) + 

(_^. s _^)2 .. c 2 (H-| K3 |F 3 (x))d 2 - c 2 9l - c 2 | K 3|(^F3(a;))a 



q,w 



0. 

(62) 
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Here we used the following notations 

L(A) = (-h 2 /2m)A + |/2| + {z/ + ^ (v) }</>o(^) + {^i + ^ (v)}#(aO> 

F 3 (s) = 2 7(1;) (<r /McQ of $(x). 

To simplify investigation of the characteristic properties of the different solutions of 
Eqs. ([5D|)-(|5TD, we subdivide the excitations fl54"D-(|55D into two major parts, the inside- 
excitations and the outside-ones. The inside-excitations are localized merely inside the 
packet area, i.e. |x| < 2L and 4> 2 (x) ~ $„, whereas the cratside-excitations propagate 
merely in the outside area, i.e. |x| > 2L and 4> 2 {x) ~ 4$^ exp(— 2\x\/L ) — > 0. 

4.1 Outside-Excitations 

For the outside collective excitations, the asymptotics of the low-lying energy spectrum 
can be found easily. Indeed, if we assume that 4> 2 {x) ~ and d x q Q (x) ~ in the outside 
packet area, the equations fl4"S| ) and ( fJ7|) are (formally) uncoupled. Then, Eq. ( p)| ) 
describes the excitonic branch of the outside-excitations with the following dispertion low 
in the co-moving frame: 

hu ex (k)^\fi\ + (h 2 /2m)k 2 , (u k (x)^n k e ikx , v k (x) « 0), (63) 

and djph(k') = q|k'| in the laboratory frame of reference. 

Then the exciton field operator, which describes the exciton condensate with one 
long-wavelength outside-excitation, has the following form: 

ip(x., t) ~ exp(— i(E g + mv 2 /2 — \fl\)t ) exp(i(y? + mvx) ) (p Q (x — vi) + 

+ exp(— i(E g +mv 2 /2— \jl\)t ) exp(i(ip+mvx) ) jexp(— i(\ji\ + his 2 /2m + k x v)t ) n k exp(ikx) > 

(64) 
It is easy to see that such a collective excitation, 

hcu cx (k) = \jl\ + h 2 'k 2 /2m + /ifc x w, 

can be interpreted as a free exciton with the energy and the (quasi) momentum 

£x(&0 = Eg + h k 2 /2m and frfcj = ft/tj + mv 5\j. 

Note that the condition huj ex > can be violated at the velocities close to v , whereas e x 
is always positive. Then the question is whether hu ex < really means the condensate 
instability in relation to the creation of outside excitations. For example, being unstable, 
the condensate could continuously emit outside-excitations, which form a sort of 'tail' 
behind the localized packet. 

Recall that the particle number N is not conserved in quantum states with a conden- 
sate, and (5N 2 ) ~ N Q . However, for N > 10 10 and T <C T c , the following estimate is 
valid 

-5 



^/(SN 2 )/N ~ 1/ y/iVo < 10" 
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Therefore, we can compare the condensate energy E (N , v) and the energy of the con- 
densate that emits excitons, or, equivalently, the condensate with outside excitations, 
(uk) ~ ySN. For simplicity's sake, we consider SN different wave vectors, {kj}, to be 
close to each other, so that the values of (k) 2 and (k x ) are well-defined. (This is a model 
of how the instability tail could be formed.) We obtain (see Eqs. fl38|)-(f4TD) 

E (N - SN, v) + E x ((k),SN) + E ph ((k'), SN) « 

+ [^-^fyN + n Cl \(k')\SN. (65) 

For the momentum of the moving condensate with the outside excitations, we have 

P ox (N Q - SN, v) + h(k x )SN + h(k' x )5N « 

« P ox (N Q , v) + (h(k x ) - mv) SN + (h(k' x ) - 3 Mv$(N , v)) SN. (66) 

Note that the energy and the momentum of the phonon part of the condensate change 
after exciton emission. We hypothesize that the transformation N —>■ N Q — 5N (with 
the emission of outside excitons, see Figs. 1,2) corresponds to the case in which the 
outside exciton and the outside acoustic phonon appear together. Indeed, in the k — > 
limit (i.e. A = 2n/k ^> L ), we approximately considered the condensate collective 
excitations as being uncoupled. However, the phonon fik' can be emitted with the energy 
compensating the changement of 8E V ^ = —(3/2) M(cf + v 2 ) i9(N , v) in the phonon part 
of the condensate energy. Moreover, the order of value of |&E7 p h| is typical for the low- 
energy acoustic phonons, ~ lmeV. If Tik' x > 0, the emitted phonon can compensate the 
changement of 5P p h x = —3Mvd(v) as well. 

The most interesting case is the backward emission of excitons, i.e. hkj = Tik 8\j < 



in the co-moving frame. Then we can rewrite (65) as follows 



E (N - SN, v) + (fuv ex (k)) 5N + (hu ph (k')) 5N w 

« E (N , v) + {2\p\ + v Q & )6N+ {\ji\ + m^f- - h\(k)\ v}6N. (67) 

The moving condensate can be considered as a stable one in relation to emission of the 



outside excitations (SN > J(5N 2 )) if such an emission gains energy, 

E (N - 5N, v) + E cxc ((k), 5N) + E ph ((k'), 5N) > E (N , v). 
This means that the following inequality has to be valid 

r h 2 (k) 2 1 

\t\ + -T^~ -h\(k)\v\ + (2|j5| + i/o € ) > 0. (68) 
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This condition can be rewritten in the dimensionless form as follows: 




2tt v 



(Am /M) Mcf/2 
\P-\(N ,v) 



>0. 



(69) 



We argue that, even for velocities close to v (where |i^(i>)| can be ~ 0.1 ^o? an d the 
instability could appear as \p,\(N , v) + h 2 kl/2m — h\k x \v < 0), inequality fl6"8|) seems 
to be always true in the long-wavelength approximation, k = i^ln/z)^ 1 < 10 -1 Lq . 
On Fig. ||, the stable ballistic condensate is shown with its long-wavelength outside- 
excitations. 

Note that the stability against large-A; modes cannot be properly described within ap- 
proximation (|63D , (|64l) . However, we can discuss this case within the inside-approximation. 



u (x 



^ (x-vt 2 ) 




1-outside-exc 



tj<t 2 v <v 



Figure 2: The ballistic condensate, 4> (x — vt) -u (x — vt)5ij, seems to be stable in relation 
to emission of the outside exciton-phonon excitations. (We consider the backward emission 
in the long- wavelength limit.) The outside-excitations presented on this figure are labeled 
by the wave vectors, k x , k' x < in the co-moving frame. In the first approximation, the 
outside-excitations can be described in terms of free excitons and free (acoustic) phonons 
emitted from the condensate coherently. 

4.2 Inside-Excitations 



To simplify the calculation of msz<ie-excitation spectrum (see Eqs. (p9|)-(pT|) ) we will 
use the semiclassical approximation ||27|| . In this approximation, the excitations can be 
labeled by the wave vector k in the co-moving frame, and the following representation 
holds: 



u 



■3,s 



x) 



a j,k 



x e 



i< / 5 k (x) 



v 3,s\ 



'i,k 



x e 



i<p k (x) 



c r a 



q k (x) 



^Vk(x) 



(70) 



where the phase </?k(x) ~ y3 + kx, and Uk(x), Vk(x), and Ck(x) are assumed to be smooth 
functions of x in the inside condensate area. Notice that the k- and x-representations 
are mixed here. This means that the operator nature of the fluctuating fields is actually 
dismissed within the semiclassical approximation. Howevever, the orthogonality relations 



among u 



3,31 



«],S'1 



and C j>s , 






and, hence, among u^k, v j,k', and C^k, 



°?,k' 



come 
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from the Bose commutation relations between the operators aij >s and oc^ s , [^fj,|27|]. For 
example, Eq. (1541) is modified as follows: 



(ty(x,t) 



dk 



u k (x)e 



i ¥'k(x)„-it<; k (x)t 



+ v^(x)e- 



-*VkW e iw k (x)t 



(71) 



and the inside-excitation part of the elementary excitation term in Eq. (|51) 
z_/j> • ■ ■ ~ 2_/j, s , out + z_/j,s,surf + J2j, s , in ' ' '■> can be written as 



Z ^ 



ls a[ on 



l.a, in 



dkdx 



(2tt) 



ftu;i k (x)ni, k (x). 



(72) 



Note that the semiclassical energy Tiu^k^) of the inside-excitation mode j, k is supposed 
to be a smooth function of x as well, (at least, as smooth as 4>l(x), which is taken constant 
in the inside-approximation). 

Although the low-lying excitations cannot be properly described within the semi- 
classical approximation, we apply it here to calculate the low-energy asymptotics of the 
spectrum. In fact, within the approximation H m H + H2, all the important properties 
of such excitations can be understood by use of the semiclassical approach. 

There are two different types of the inside-excitations, the longitudinal excitations 
and the transverse ones. The later have the wave vectors k perpendicular to the x (in- 
direction. For the sake of simplicity, we choose k || Oy. Then the vector CJ' k has one 
nonzeroth component for such transverse excitations, C v , k 7^ 0. 

Substituting ansatz flTOD with k r = k±52, r and Cj k = C V j k b^ r into Eqs. (p5[)-(|5T]), we 
transform these differential equations into the algebraic ones (within the inside-approximation 
L(A)^L(-k 2 )): 



(L(-k x ) - %0J i>k ) u fc (x) + \u <f> o (x) + 2u 1 (j) (x)j Vj- ifc (x) + (r (j) o (x) ifc ± C| fc (x) = 0, (73) 
uo4>l{x) + 2i/i$(z)) u i)fc (x) + (L(-kl) + hu j:k ) Vj-, fc (x) + aoM%) ik±CK ik (x) = (74) 



P vo<fio(x) ik± Uj,fc(x) + p a (j) o (x) ik± v i)fc (x) + 



u> 



2;„2 



3,k 



Cfk 



c-> 



j,k\ 



0. 



(75) 



After some straightforward algebra, we can write out the equation that defines the spec- 
trum of transverse exciton-phonon excitations in the inside-approximation: 



(ftu; 



3,k) 



(ujl k -c*k ± 2 )x 
L(-ki) - u <f> 2 o (x) - 2u 1 <f>i{x) ) (L(-k 2 ± ) + uotKx) + 2u x <j>^x) 
L(-k 2 ± ) - u^Kx) - 2u 1 <f>i(x) 



2 -4<Pl(x)(ctkl). 
pcf 

Taking into account the momentum cut-off ko, which is defined as 

(h 2 /2m) k 2 « \jl\ = (h 2 /2m) L^ 2 , k > k , 



(76) 
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we can rewrite Eq. Q7B|) as follows 

("Ik ~ tf k ± 



X 



X 



(hu^) 2 - (^-[kl - kl] + F{x) + e\ ■ (^-[kl - Q + F{x) + 2vufi(x) + e' + 



= l—[kl-k4]+F(x) + e + 2j^(a at)$i(cikl), \x\ < L . (77) 

Here F(x) = \vq(v)\(Q>q — <fil(x) ) > 0, i.e. F(x) — inside the condensate, and 

e + ^z^0^(x), e' + « 5i/i0„(x), k± > k . 



Although Eq. ( [77]) can be solved exactly for the transverse excitation spectrum [30], 
taking into account the coupling term in the r.h.s. of flT7| ) changes the values of excitation 
energies slightly, and the excitations can be approximately considered as of the pure 
excitonic {Tiuji^ = hu eXjk± ) or the pure phonon {Ti^k = ^ph,k±) types. 

It is also useful to investigate asymptotics of the transverse inside excitations. For 
defmiteness sake, we investigate the left side asymptotics of these excitations here, 

u i)S (x) = exp(l a x/L ) e ik±y u j>k , v iiS (x) = exp(e v x/L ) e ik±y v j>k , x<0, (78) 

q s (x)=exp(4x/L )e ifc ^q fc , q,(x) = exp(4x/L ) e ik ^ y q fc , x < 0. (79) 

Here Uj t k, Vj,k, Cj^, and C* k are smooth functions of x at \x\ > L . Note that we 
introduced two components of C!j s ~ e lk±y to make Eqs. (|59|) -(|62"D self-consistent. Let the 
equalities 

£ u = £ v = l and 1 + £ u = £ c = 2 (80) 

be valid. Then the system of differential equations fl59"|)-(|B2"D can be reduced to a system 
of algebraic ones, which are analogous to Eqs. (|73|)-(|75D. Consequently, we can write out 
the equation for u)j tk {x) valid at |x| > L , 

( [u jtk - iv (2/L ) f - c]\k\ - (2/L ) 2 ] ) x 

(H>) 2 - (£^[kl - *g] + H*)) ■ (^l ~ ~ k l\ + H*) + 2^o(2$ exp(x/L ) ) 2 ^ 

= ^[kl-~k 2 }+F(x)^j 2^(a af)(2$ exp(a;/Lo)) 2 c 2 [A;i-(2/L ) 2 ], (81) 

where fco = V%/Lo — k , k± > k , and 

F(x) = \Vq(v)\ ($ 2 - (2$ o exp(x/L )) 2 ) -» 2|//| -e at \x\ > L . 

(We neglected the terms, such as v\ <\)\{x) % )S (x) and v\^ [x) Vj iS (x) ~ exp( (4 + £ u )x/L ) 
in Eqs. (|59|) -(|60D, and the terms oc k 3 in Eq. fl62p as well.) 
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Obviously, the structure of Eqs. ( [77| ) and (^) is the same. As the coupling between 
exciton and phonon branches is weak for the transverse inside-excitations (see the r.h. 
sides of Eqs. (|77|) and fl8lD ) and the effect of the finite width L can be taken into account 
as the spatial dependence of the important parameters in hu ex , we use the following 
formula to estimate the low-energy excitation spectrum: 

(hu ex>k± ) 2 ~ (Jj^(kl - kl) + F(x) + e + \ {^Akl - kl) + F(x) + 2v <f>l(x) + e' + ) ~ 

h 2 

~ — (k 2 x - kl) 2z/ $2 + 2z/ $2 e+ at k x -» fco. (82) 

2m 

Here we take the inside-condensate-asymptotics of F(x), (fil(x), and e + ~ v{(v) 4>t( x ) to 
estimate /k^ ex . Note that, for the inside-condensate excitations, the the low-energy limit 
means 

(h 2 /2m)(k\-kl) ~(l~10)|/i|. 

Then, in the co-moving frame, the low-energy excitation spectrum hu ex k ± niay develop 
a gap of the order of \p,\, (see Eq. ( |64[) for comparison). Thus, inside the condensate, 
we obtain a strong deviation of the collective excitation spectrum from both the simple 
excitonic one, |//| + (h 2 /2m)k 2 _, and the Bogoliubov-Landau spectrum oc \k±\. 

4.3 Longitudinal inside-excitations 

The case of the longitudinal excitations, k r = k x 6i >r , Cj k = C? k 8i >r , is more difficult 
to analyze because the mode interaction is non- negligible in the low-energy limit. (On 
Fig. 1., a longitudinal inside-excitation is shown with the two possible directions of the 
wave vector k || Ox.) Recall that the "bare" phonon modes, which can be written in the 
laboratory frame as 

u x (x, t) ~ q (x — vt) + Cl(x — vt) exp(ik x x — iujp^t) + c.c. 

with Wph = ci\k x \ and Cf(x) ~ (fil(x), will be considered in the co-moving frame, x — vt — ► 
x. Then, within the inside-condensate approximation, the following equation stands for 
the excitation spectrum: 

( (u jt k + vk x ) 2 - c 2 k 2 x ) x 

(^, k ) 2 - (^[k 2 x - k 2 ] + F(x) + e + \ ■ (^[k 2 x - k 2 ] + F(x) + 2M 2 {x) + e'+ 



g-[k 2 x - k 2 ] + F(x) + e + j 2j^ (a af) ^ (c 2 k 2 x ), \x\ < L . (83) 

It is important to note that, unlike the case of transverse excitations, the values of 

(fuv^ kx f^h 2 ( Cl -v) 2 ((3^7)k ) 2 



22 



and 

WW" * (|^( 10 - 4 °) fc o + * + ) (^(10 ~ 40)A; 2 + 2u 0( pl(x) + e' + ) 

are of the same order of value at k x ~ (3 ~ 8) fco, and the inequality (hu^ kx ) > (huj p ^ kx ) 
is valid in the low-energy limit. Moreover, the two cases, k x > (+-case) and k x < 
( — case), are different as it can be seen from the l.h.s. of Eq. (|83|). In the low-energy 
limit, we can write the approximate solution of (J83| ) as follows 

041) 2 « (^[^ - fc o] + F{x) + e + ) x 
x f ^{k 2 x - kl] + F(x) + 2u <f)l(x) ± 2 g ±7 (v) ^ (a af) <f> 2 Q (x) + e' + J , 

,(+) ^ fi,.,(0) 



where q + ~ 1 and < g_ < 1. Note that hu eyik > fiuJ c J k , whereas, for the phonon-type 
,(+) 

J ph,l 



branch, aA k < (q — v ) fc x . For k x < 0, we have the following inequaluty for the excitonic 



branch, 

(h 2 



{) Jkl-k 2 }+F(x) + e + \x 

x {^[k 2 - k 2 ] + F(x) + 2V,{v) <f> 2 (x) + A < (Mx",l) 2 < H? fc J 2 , (84) 

where 2uq(v) <pl(x) ~ — 4|/2| within the inside-approximation, and, for the phonon-type 
branch, we obtain oj p ^ kx > (q + v ) \k x \. 

To derive the formulas for the amplitudes u k (x), v fc (x), and C k (x) of the excitonic 
branch, we use the following approximations 

M-© = L(- k D + ^-gM^ M L( _^ ± g±7(u) _^ (wt?) ^,, 

and B = v Q (f) 2 (x) + 2z/i0q(x) is modified as 



B± « z/ ^(x) + 2i/i0*(x) ± g±7(v) -77^ (a af ) #0*0 ■ 



Then we can rewrite the formulas for the excitonic excitation spectrum as follows 
(&"Sj a « 4H£) - B 2 ± = (L(-k 2 x ) - B) (L ± (-k 2 x ) + B ± ). 



Recall that the orthogonality relation (|5q) can be used to normalize the amplitudes. 
Within the inside-approximation, Eq. (E^) can be rewritten as follows (5 SS = 1 — ► 5^ = 1) 

dx(|u fc (x)| 2 -|v fc (x)| 2 ) + (l/K)jd*2p{u^ kx +vk x )\Ct(x)\ 2 = l, (85) 
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and we have for the excitonic amplitudes 

2 /T,\ T^(-k 2 )-huj {± \ 

(86) 



VL^ix) 



2 fr±\L±{-K )+hu eXtk 



-' * .' = ' 2 /t ± u ± K)-^» 



r^_j 






Here the effective condensate volume V^s ~ 2S"Lo is used to normalize the u- and v-wave 
functions of the inside excitations, and fdr (\n k \ 2 — |v fc | 2 ) = T± < 1. 
Subsiquently, we get for 

Qxt x s = P'Vq <Pq(x) ik x (ufc(ar) + v k {x) ) 

" (uj cx , kx + vk x ) 2 -cfk 2 

the following approximate formulas 

C x k {± \x) « T g ±7 (^) ^\A^oWf yfi${uF>{x) + vi ±} (x) ). 



To estimate the characteristic value of C^ (x), we use Ufc(x) — v&(a;) ~ \J"£±/V c k and 
obtain 

yx(±)i _. _/_\ a o /. s ^, lY±aiL ai 



The parameters T± characterize the relative weight of excitonic degrees of freedom in 
the considered branch of excitations. As the parameter hv k x /fiuj ey . kx < 1 at k x cz (4 ~ 
8) ko, the parameter T + (k x ) can be estimated as 0.5 ~ 0.7. For k x < 0, we obtain the 
following equation from (|85*|). 

/ (qf «&g) , , , crp ,2 ( %y\k x \ \ L{-kl) - B 
[ + m/M 9 - 7W Wc?J ^ fo^J ft a *g/2ro 

Within the stability area (see the next subsection for an extended discussion), we estimate 
the ratio hv\k x \/hw^ kx ~ 1/2 ~ 1/3 at (^1 ~ (4 ~ 8) k . Then, T-(k x ) > and 
T_~0.6~0.8. 

To go beyond the inside-approximation, the effect of inhomogeneous behavior of the 
longitudinal excitations can be considered. We use the following ansatz for the left side 
asymptotics (see Fig. 2) 

u i>s (x) = exp(4x/L ) e ik * x u k , v 3> (x) = exp(£ v x/L ) e ik * x v k , (88) 

q s (x) = exp(4a;/Lo) e lkxX C%, x < 0. (89) 

where £ u = £ v = 1 and £ c = 2. Then, like the case of transversal excitations (see Eq. 



81)), we can write the equation for Uj :kx (x) valid at \x\ > L , 

((^j,k x +vk x ) 2 ~ cfk 2 x ) x 
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X 



(S^fc,) 2 - (|^K - *o 2 l + m) (|^E - *§] + F(x) + 2„„(2* o exp(x/L ) ) ! 



^k-^m)^. 



2 ;„2l , Bf„\ o °0 C_ „3\ /o* ™„^„/r \\2„2k2 



[A-: - A;;] ./•'(,•) j 2^ (a af) (2$ exp(x/L ) Y cfk^ (90) 



where k x ^ k x = k x — i(2/Lo) in the phonon parts of this equation and k x — > fc a = 
fc x — i(l/L ) in the exciton parts of it (x < 0). It is easy to see that Eqs. ( JS3D and (|90|) are 
in the continuity correspondence, i.e. they describe the same object. For example, the 
(left side) asymptotic behavior of Uoj^ k (x) can be obtained from the inside-condensate 
formulas by the substitute 

k x ^k x , F{x) -► F(x) -► 2|/2|, 

and 

02( x )^(2$ o exp(x/L o )) 2 ^O. 

As a result, we obtain fuv eXt k x — \p>\ + ti 2 k 2 c /2m that corresponds to the outside-excitation 
spectrum, Eq. (|63|) . 

The excitonic input into (5$^ 5$ (x))t=o, the quantum depletion of the moving conden- 
sate, can be calculated by J2i, s |vi, a (a;)| 2 p7| - To estimate this value, one can approximate 
|vfc(x)| 2 as follows 



v£ } (x) 



2 /T ± \ 5 2 



r^_/ 



,Kff/4L 2 ± (-^)- 

However, the summation J2i s implies Jdk x d 2 k±/(2ir) 3 within the semiclassical approx- 
imation. Assuming that such an integration makes the difference among vjf" , vjj. and 
v fcx not essentially important, we use the following estimate for |v k (x)| 2 , 



v k(£, 



2 / 1 \ B 



r^j 



i 



K ff /4L 2 (-k 2 ) 



Then, the integration /<ifc x d 2 k±/(2n) 3 can be reduced to J k k 2 dk/27T 2 , and the main input 



(~ <Po( x )) can be estimated from the following formula 

/X&W( \\ ~ l W&Q i (Ml(x) f 1 vo$ 2 cosh- 2 (x/L ) 

{ l ))T=0 "8vr 2 L3 2|/2| 8vr 2 L3 2|/i| 2 ~ 8tt 2 Lg N $ 2 

Using this density, we speculate that the localized depletion of condensate - i.e. the 
number of particles that are out of the condensate but move with it coherently - seems to 
be a small value. We obtain the following estimate (z/ = s d 3 B and |^o(f)| = l^oMI &%)'■ 

6N = [dx(6&6$tx))T=o * t^AtI^ 

J lu IT // J-jq 

where the factor before 7V can be estimated as 



"h . *oMv)\ 1 (10 -i^ 10 - 2) - : 



\p,\(N ,v) L 3 (N o ,v) 2x(Ry*) 2 n v ' ° 
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(Here we used Eqs. Q55D , (|5fi|) within the approximation n 3> 10.) Note that there is 
no small k input to the estimate of J2i, s because, first, such excitations belong to the 
outside-excitation branch in our model, and, second, we use the approximation Vk(x) ~ 
for them. 

4.4 Stability of the Moving Condensate 

To investigate the stability of the moving condensate in relation to the creation of inside- 
excitations, we calculate the energy of the condensate with the one inside excitation, 
(a\ s ai tS ) = 1, described by the following set: k, u k , and u&, v fc , C k . In this study, we 
analyze the stability inside the excitonic sector of our model. 

Although the excitations were defined in the co-moving frame, calculations should be 
done in the laboratory frame. Returning to the lab frame, we represent the exciton and 
phonon field functions as follows: 

4> (x — vt, t) — > (f> (x— vt, t)+exp(— i(E g + mv 2 /2 — \jl\)t ) exp(i((p+mvx) ) 5\&(x, t), (91) 

where 

<5§(x, t) = U k (x - Vt) e*^ +kx ) e -^ k +k*v)t + Vfe ( x _ yQ e -*(^o+kx) e i(u, k +k x v)t^ 

and 

u (x — vt) —* u (x — vt) + Ck(x — vt) exp(i(<£> + kx)) exp(— i{u k + k x v) t) + c.c, (92) 

see Eq. (^) for comparison. In this analysis, the inside excitations are not considered as 
fluctuations, and the (average) number of particles in the condensate and its energy are 
changed as N Q — Jdx.5ip^5ip and E (N ) — (<9at-E ) fdyLd~ip^5ip, respectively. However, these 
changes are not important if the number of excitation in a system is less than VSN 2 ~ 
y/Np. They could be important in the case of instability of the moving condensate. 

The zeroth component of the energy-momentum tensor can be represented in the form 

T ° = T O °(0 O , Uo ) + r 0(2) (^t ; s*, Su, d t 5u | O , Uo ), 

where the first part corresponds to the condensate energy E and the second part gives 

r o- 

'-o 



the energy of inside-excitations, E- m . After substitution of (|91p,(|92]) into E in = fdx.T 



we have for the total energy 

J - J o T^ -'-'in— ex 

« E (N ) + 5E (N ) + Jdxh(uj k (x) + k x v) {\u k \ 2 - |v fc | 2 + (2/%)p(uj k (x) + vk x )\C k \ 2 ) = 

= E (N ) + 5E (N ) + h(u k + k x v), (93) 

where 

6E (N ) = (2|/2| + v $ 2 o ) [dx6ft6i> - M ^ +v ^ 3 tf(iV , v) [dx6ft6#, (94) 
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SftS^x) -> |u fc | 2 + \vk\ 



2 



see Eqs. (|6q),(p7|) for comparison. 

In this study, we discuss qualitatively the stability of the condensate in relation to the 
backward emission of inside-excitations, (i.e. k x < in the co- moving frame). To begin 
with, we consider the standard criterion, 

H u ix,k — \k x \v)>0 at \k x \~ zLq 1 , (95) 

where z ~ 3 ~ 10 corresponds to the low- lying excitations. The value of lu cx k (x) is taken 

within the inside-approximation, see Eq. (|8~4f) , so that huj cx k ~ f(z) \ji\ and z/f(z) ~ 
0.1 ~ 0.3. Then, it is easy to conclude that the following inequality 



h\k x \v 



(4m /M) Mcf/2 , s 

< 1 (96) 



is valid in the low-energy limit if the effective chemical potential Q3TSD is large enough, see 
Eq. (|69"D for comparison. 

More precisely, the ballistic velocity v and the number of particles in the condensate, 
N Q , have to be large enough, f.ex., |£o(u)| — (10 _1 ~ 1) Ry* and n ~ 10, in order to the 
inequality 

{4m/ M) Mcf/2 _ (4m/M) Mcf/2 

\fi\(N ,v) ~ \s (v)\y4nlxRy* ~ l } 

can be satisfied. Thus, for 

\fx\(N , v) > /i cr ~ 10- 1 {4m/ M) Mcf/2, (98) 

where // cr ~ 10 _4 eV, one can expect conditions (p6|),(p7|) to be valid. 

Despite the condensate can be formed near 7(f) w 7 in theory, f.ex., with 1^(^)1 — 
0.1 Vq and n Q ^> 10, such a ballistic state seems to be unstable against the creation 
of inside-excitations. Note that the critical (Landau) velocity, t> cr , can be found as a 
solution of Eq. ( p8|) and v < v CT (N ) < c\. In fact, the parameter \jl\/ ji cr controls the 
stability/instability of the condensate, see Eqs. (|69"D , (|96|) . 

Analyzing (|95]), we did not take into account Eq. (0). However, if the instability 
regime takes place, more than \fNl inside-excitations can appear. As the changes in u (x— 
vt) because of N — > iV — 5N are nonlocal (in spite of creation of the localized excitations, 
see Figs. 1 and 3), a free acoustic phonon can appear in the system lattice + excitons 
together with appearance of the localized excitation huj c J k {x). Like the case of outside- 
excitations, we assume that (3/2) M(cf + v 2 ) $(iV , v) ~ hcik p h, see Eq. fl9"3f). Then, only 
the term oc 2|/t| + vo&f, is important. In fact, this term leads to some renormalization of 
the values of the critical parameters, /x cr and v CT . 
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Figure 3: The ballistic condensate, (f) (x — vt) ■ u (x — vt)5ij, can be unstable in relation 
to emission of the inside-excitations if the effective chemical potential \ji\(N , v) < fi cr . 
In terms of Landau critical velocity, this means v < v < v CT (N ). If such an instability 
takes place, the emission of inside-excitations can be accompanied by the emission of 
outside-excitations of the condensate. The longitudinal inside-excitations are labeled by 
the wave vector k x < on this figure, whereas the outside-excitation is labeled by the 
wave vector k' x < 0. 



5 Interference Between Two Moving Packets 

In this section, we address the problem of interaction between two moving condensates. 
This problem is essentially nonstationary, especially if the initial ballistic velocities of 
packets are different. Within the quasi- ID conserving model, the following equations 
govern the dynamics of the two input packets (we choose the reference frame moving with 
the slow packet): 



,^2 



ihd t H — )VoOM) = (—z—dl + vo\ipo\ 2 + vi\ipo\ 4 )*Po(x,t) + a d x uo(x,t)i/j (x,t), 



(99) 
((d t -v'd x ) 2 -cfdl)u (x,t) -cf 2K 3 <%uo(x,t) d x u (x,t) = p~ l a d x \i> \ 2 (x,t). (100) 
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Figure 4: Two ballistic condensates move with different velocities, u — 1>' ~ (O.f ~ 0.3) q, 
and t — t\ before the "collision" , or, the strong interaction process. If one can prescribe 
the coherent phase to each of the participating condensates, e.g., <£> c (x) « tp + (mv/h) x, 
the interference area appears between them. (The interference area is marked by bold 
dashed lines on this figure.) As v ^ v', the fringes are non-stationary, and the outside- 
excitations can actually be excited in this area. 

Then, the initial conditions can be written in the explicit ID form by using the exact 
solution of the model (PH|) , (|Tl~D . Note that the amplitudes of the stationary ballistic state, 
(j) (x — vt) -d x u (x — vt), were defined from the normalization condition and depend on the 
values of v and N . Hence, the amplitudes of the "input" condensates for Eqs. (|99f), (|100|) 
may not have the same values. 

In this study, we approach the problem of strong interaction between the condensates. 
Therefore, we choose the nonsymmetric initial conditions, i.e. the amplitude and the 
velocity of the "input" packets are different, for example, v > v' and v||v', see Fig. 
f|. Here, we reply on the experimental observation || that, at v > v , the ballistic 
velocity of the condensate depends on the power of a laser beam irradiating the crystal. 
If the exciton concentration in the first packet, n \, is close to the value of the Bose 
condensation threshold, the exciton concentration in the second packet, n o2 > n ol , the 
velocity difference between condensates can reach (0.1 ~ 0.3) q. Then, in the reference 
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frame moving with the first (slow) packet, the initial conditions can be taken as the 
following: 

if>o(x,t = 0) • u (x,t = 0) = 0oO; -/Vol) ■q (x;N ol ) + 

+ exp(i(5(p + mdvx)) (f) (x + x ; N o2 ) ■ q (x + x ; N o2 ), (101) 

where S(p — tp — tp', Sv — v — v', xq = v't, and r is the (initial) time delay. As the 
second packet moves in this frame of reference, the regime of strong nonlinear interaction 
between the condensates is (theoretically) unavoidable. Note that, even before collision, 
a time- dependent interference term in \tp (x,t)\ 2 begins to influence the packet dynamics, 
see Fig. 4. For example, the r.h.s. of (|100|) contains 

~ d x { 2 cos(m 5vx — 5ut + 5(p) 4> Q (x; iV ol ) (fi (x — 5vt + x ; N o2 ) } , (102) 

where 

h5uj = m5v(v + v')/2-(\fi\(N o2 ,v)-\fi\(N ol ,v')) 

and |/t| oc ( N Q /N*) 2 \io(v)\ 2 . The ratio | /5-2 1 / 1 /xi | can be of the order of 10 1 , and the 



characteristic scale of fringes ( 102 ) is 
7r(h/m5v) 



(10 ~ 30) 



lflW.i,«> ,. Bm10i 



L V J \ (m/M)Mcf/2 

that is they are of the long-wavelength nature. 

To answer the question which model (conserving (|99|) , ( |100] ) or kinetic £L9|.|32|] one) 



is more adequate to describe the packet collision, we have to compare the estimate of 
interaction time, 

r* =* L ch2 /Sv ~ 10 3 a B /0.2q ~ 10" 9 ~ 10" 10 s, 

and characteristic time scales of the processes 

|iVo2(i),^>-|iV i(£)y> - \N o2 (t)±5N,v)-\N ol (t)T5N,v') (103) 

driven by phonons or by x — x interaction. Note that some thermal phonons have to be 
excited in the system to assist such transitions, and the value of r* is of the order of scat- 
tering time of the exciton-LA-phonon interaction (although without any macroscopical 
occupancy) |L7|. 

If processes ( |103j ) are driven by the lattice phonons, two phonons are necessary to 
satisfy the laws of conservation. For instance, we choose SN = +1 in (|103|) and obtain 
(see Eqs. 



hki, x = m5v + 3Mvtf( y N o2 , v) - 3Mv' tf(N ol , v') + Tik 2 ^ x , 

hctlhj = m5v (v + v')/2 - 3 ( |/2|(iV o2 , v) - |/2|(iV ol , v') + v Q $^ 2 /3 - u $ 2 ol /3) + 
+ 3/2 M{c 2 + v 2 ) $(N o2 , v) - 3/2 M{c 2 + v' 2 ) &(N ol , v') + % Cl \k 2 J. 
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Although the second packet moves faster, mSv > 0, this state can be considered as a more 
stable (and, thus, more preferable) one for the excitons of the slow packet. Indeed, the 
following inequality for the effective difference between the generalized chemical potentials 
seems to be valid 

mSv (v + v')/2 - 3 ( \jl\(N o2 , v) - \jl\(N oU v') + u $* 2 /3 - u & ol /3) < 0, (104) 

and the absolute value of the l.h.s. of ( |104p is ~ \p,\(N o2 v). Thus, within the quantum 
kinetic model, the relevant transition probabilities have to be calculated at least in the 
second order of perturbation theory, f.ex., 

\N o2 ,v)-\N ol ,v'} -^ \N o2 ,v)-\N oU huj cx , kj ,v f } ^ \N o2 + 2,v)-\N ol -2,v'). 

As a result, the system of two quantum Boltzmann equations will describe the interaction 
process. 

Unlike the Boltzmann equations, Eqs. (|99| , |100| ) contain information about the quan- 
tum coherence between two condensates explicitly and, moreover, can describe the case 
of strong interaction. Unlike the ID NLS equation that supports many-soliton solutions, 
Eqs. (|9^),( |100| ) are, in fact, quasi-lD ones, and u > in (^). Therefore, it is an open 
question what happens with two (excitonic) solitons after they collide in the crystal. 

In this work, we assume that the dominant process(es) of condensate interactions 
is that one(s) leading to d t N o2 > 0. Then, at the time scales ^> r*, one solitonic 
packet can appear as a result of these processes. Such a resultant ballistic packet can 
be approximately described by the steady-state one-soliton solution of Eqs. ©,(11) with 
N Q ?« jV 2 + N 1 and the low of energy conservation, 

E(N ol ,v') + E(N o2: v)nE(N = N o2 + N ol: v). (105) 

If N < N*, all the approximate solutions having been found in this study are valid to 
describe the resultant packet. 

As we prescribed the value of N Q , we have to estimate the value of v from Eq. ( |105| ), 



(generally, v ^ v ). Moreover, we have to assume that the total momentum of the conden- 
sates, P x (N i,v') + P x (N o2 ,v), may not be conserved because of lattice participation in 
such a condensate "merger" . However, the challenging question of the results of coherent 
packet collision needs further theoretical and experimental efforts. 

6 Conclusion 

In this study, we considered a model within which the inhomogeneous excitonic condensate 
with a nonzero momentum can be investigated. The important physics we include in our 
model is the exciton-phonon interaction and the appearance of a coherent part of the 
crystal displacement field, which renormalizes the x-x interaction vertices. Then, the 
condensate wave function and its energy can be calculated exactly in the simplest quasi- 
lD model, and the solution is a sort of Davydov's soliton |23|. We believe that the 
transport and other unusual properties of the coherent para-exciton packets in Cu 2 can 
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be described in the framework of the proposed model properly generalized to meet more 
realistic conditions. 

We showed that there are two critical velocities in the theory, namely, v and v CT . The 
first one, v , comes from the renormalization of two particle exciton-exciton interaction 
due to phonons, and the bright soliton state can be formed if v > v Q . Then, the important 
parameter, which controls the shape and the characteristic width of the condensate wave 
function, is |/i|/V\ see Eqs. (^),(p6|). The second velocity, t> cr , comes from use of 
Landau arguments |2B| for investigation of the dynamic stability / instability of the moving 
condensate. In fact, the important parameter, which controls the emission of excitations, 
is |/i|//i cr , see Eq. (|98|) . Then, within the semiclassical approximation for the condensate 
excitations, we found more close v is to q (/x cr < \jl\(N , v) < /i*) more stable the 
coherent packet is. It is interesting to discuss the possibility of observation of an instability 
when the condensate can be formed in the inhomogeneous state with v 7^ 0, but with 
v o < v < v CT (N ), or, better, \jl\(N , v) < /i cr . Such a coherent packet has to disappear 
during its move through a single pure crystal used for experiments. As the shape of the 
moving packet depends on time, the form of the registered signal may depend on the 
crystal length changing from the solitonic to the standard diffusion density profile. 

We found that the excited states of the moving exciton-phonon condensate can be 
described by use of the language of elementary excitations. Although the possibility 
of their direct observation is an unclear question itself (and, thus, the question about 
the gap in the excitation spectrum is still an open one), the stability conditions of the 
moving condensate can be derived from the low-energy asymptotics of the excitation 
spectra at T <C T c . However, the stability problem is not without difficulties [^0|] ,|33|. 
One can easily imagine the situation when the condensate moves in a very high quality 
crystal, but with some impurity region prepared, f.ex., in the middle of the sample. 
In this case, the excitonic superfluidity can be examined by impurity scattering of the 
ballistic condensate. Indeed, such impurities could bound the noncondensed excitons, 
which always accompany the condensate, and could mediate, for instance, the emission 
of the outside excitations. The last process may lead to depletion of the condensate and, 
perhaps, some other observable effects, such as damping, bound exciton PL, etc.. 
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